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Abstract 

A continuum model for a population of self-propelled particles interacting through 
nematic alignment is derived from an individual-based model. The methodology 
consists of introducing a hydrodynamic scaling of the corresponding mean-field ki¬ 
netic equation. The resulting perturbation problem is solved thanks to the concept 
of generalized collision invariants. It yields a hyperbolic but non-conservative sys¬ 
tem of equations for the nematic mean direction of the flow and the densities of 
particles flowing parallel or anti-parallel to this mean direction. Diffusive terms 
are introduced under a weakly non-local interaction assumption and the diffusion 
coefficient is proven to be positive. An application to the modeling of myxobacteria 
is outlined. 
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1 Introduction 

Systems of large numbers of self-propelled particles are commonly found in nature. In 
three space dimensions, well-known examples include hsh schools, flocks of birds and in¬ 
sects [Hl[inil39]. In two space dimensions, typical examples are bacteria, such as myxobac¬ 
teria and Bacillus subtilis [SlESlilliniEH]. However, the interest in such systems is not 
limited to biological examples. It has been found that the underlying principles also often 
apply to human walking behavior [HUl ES] and physical systems, such as granular media 
[2B1 IS7j . The main interest in these systems, both mathematically and from an applica¬ 
tion’s point of view, stems from their emergent properties. Emergence refers to the fact 
that fairly local interaction rules between the agents like attraction, alignment, repulsion, 
etc. give rise to large scale collective behavior. Emergence is physically described as a 
phase transition to an ordered state (see e.g. the review [H]). Mathematical models often 
rely on a particle (agent) based description [H, HU [I31I1I] • Among these models, the time 
dependent Vicsek model m and variants of it have received much attention owing to its 
structural simplicity, while convincingly describing a wide range of phenomena. The gen¬ 
eral idea of the Vicsek model is to describe the agents as particles, moving with constant 
speed. At each time step each particle adjusts its direction towards the mean direction of 
the particles in its neighborhood subject to some random perturbation. 

Whilst particle models are a good tool to reproduce and study the observed phe¬ 
nomena, the high number of agents involved make simulations computationally costly 
considering the amount of detail one is interested in. Hydrodynamic (macroscopic) mod¬ 
els are typically easier to simulate and analyze and thereby offer a powerful alternative. 
Therefore a lot of effort has been made to formulate macroscopic equations starting from 
particle models P El [9l [121 EU ESI EEl IBS]- Mathematically one of the main challenges is 
the lack of conserved quantities, which is typically a key ingredient in the computation of 
the macroscopic limit (see again [15]). This lack of conservation is frequent in biological 
systems, which, as opposed to many physical systems, e.g. don’t conserve momentum. 
One strategy to overcome this is the use of generalized collision invariants (GCI), intro¬ 
duced for the hrst time in |23] to derive a hydrodynamic model of the Vicsek dynamics 
(later on referred to as the Self-Organized Hydrodynamic (SOH) model). The passage 
from individual-based to hydrodynamic models involves the construction of an interme¬ 
diate kinetic model. For classical physical systems, momentum or energy conservation 
translates at the kinetic level into the concept of collision invariants. This is general¬ 
ized to the GCI concept if the underlying system does not satisfy these conditions. The 
rigorous passage from kinetic to SOH models has been achieved in |36]. The GCI con¬ 
cept has been used in a variety of cases HU 11611201128] and even applied to a model of 
economics [22] • Modal analysis of the SOH model can be found in [23]. One important 
feature of the Vicsek model is the emergence of phase transitions. These are studied in the 
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kinetic of macroscopic framework in [21 [6], [181 [19], US] . Here, we concentrate on the deriva¬ 
tion of macroscopic models and postpone the analysis of phase transitions to future works. 

In this paper we want to apply this method to nematic alignment of polar, rod-shaped 
particles moving with constant speed in two space dimensions. Nematic alignment de¬ 
scribes the situation where upon meeting, two particles either align if their mutual angle 
is acute or anti-align if their angle is obtuse. This nematic alignment rule typically results 
from volume exclusion and can be understood as an inelastic collision between polar par¬ 
ticles. It differs from polar alignment of polar particles (which is akin to ferromagnetic 
interactions in spin systems, as is the case of the Vicsek model) and nematic alignment 
of apolar particles (active nematics, which is typically observed in liquid crystals [H], 
in hbrous tissues [T3] or in systems of interacting disks [21] )• Large collections of polar 
particles which nematically align have been observed to exhibit interesting macroscopic 
patterns such as high density traveling bands and their properties differ from the other 
types of alignment [29l [38] . Therefore the need for a hydrodynamic model derived from 
particle dynamics arises and this work suggests how to derive such a model. In the deriva¬ 
tion we account for purely local interactions as well as weakly non-local interactions which 
leads to an additional diffusion term in the macroscopic limit. 

To give an example how the nematic SOH model can be applied to biological systems, 
we consider the case of myxobacteria. Due to their ability to produce several macroscopic 
structures, such as waves, spirals and clusters [271 [33l [16], they have already received a lot 
of attention in terms of mathematical models [3 EB 1321 ST]. Since this type of bacteria 
can internally reverse its polarity, the model has to be complemented with a reaction 
term, describing these reversals. We demonstrate how such a term can be incorporated 
in a systematic way into the original model. 

The rest of the document is structured as follows: In Section we describe the particle 
model for nematic alignment and the resulting mean-held model. The derived alignment 
operator is examined in Section Whilst it shares some properties with its polar 
alignment counterpart, there are also some differences, e.g. the set of equilibria has 
dimension three (instead of two). Mathematically more unexpected, showing the existence 
and uniqueness of solutions of the related linear operator, which amounts to proving a 
Poincare inequality, is slightly more complex than in the polar case [23] • The section 
hnishes with explicitly characterizing the equilibria and generalized collision invariants 
of Q^i- With these ingredients we derive the macroscopic model in Section]^ and prove 
its hyperbolicity. In Section we assume weakly non-Iocal interactions and derive the 
corresponding macroscopic model. Finally in Sectionj^we apply both macroscopic models 
to the case of myxobacteria by additionally describing density dependent reversal at the 
particle level. In the macroscopic model this leads to a reaction term which structurally 
hints at interesting dynamics expected to take place. 
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2 Particle and Mean Field Model 


2.1 Particle Model 

As a first step towards a macroscopic model, we start with the description of a time- 
continuous particle (or agent-based) model. The model’s structure is similar to the time- 
continuous version of the Vicsek model presented in [23] and we will discuss similarities 
and differences. We image the particles to resemble stiff rods of equal length and describe 
the Ath particle {i G A/" := {!,...,iV}) at time f > 0 by its center of mass Xi{t) G 
and its orientation n(0j(f)), where we define 

\ smO J \ cosO J 

as the polar vectors associated with the angle 0 in a fixed reference frame and where 
0j(f) is defined modulo 27r. All particles are assumed to move in the direction of their 
orientation with constant speed, Vq- In particular this means that they cannot drift side¬ 
ways (as opposed to the model in [1]). Under these hypotheses, the motion of the particles 
is supposed to satisfy the following system of stochastic differential equations 

d X 

^=Vov(e,(t)), ( 2 . 1 ) 

d0j = —z/Sign( cos(0j — 0*)) sin(0j — 0j) dt + (2d cos^(0j — 0j))^'^^ di?^*. (2.2) 

Here 0j = Qi{t) is the average nematic alignment direction. It is defined as follows: First 
dehne a mean ’nematic’ current 


M't) = v(20fc(f)), 

keJV:\Xkit)-Xi{t)\<R 


where i? > 0 is the interaction radius. Next dehne Qi{t) as an angle of lines (i.e. an angle 
dehned modulo tt) such that: 


!i(2ei(«)) 


Mt) 

Mt)\' 


(2.3) 


We note that the dehnition of the vector on the left-hand side of (2.3) does not depend on 
the choice of the reference frame used to measure the angles 0* as Qi{t) can be equivalently 
dehned by 

^ sin2(0fc(t) - Qi{t)) = 0, 

keJV:\Xkit)-Xi{t)\<R 


and we notice that this relation is invariant under a uniform translation of all angles. To 
obtain a unique representation of the angles, we will use the convention that the particle 
angles Qi{t) are taken in [— tt, vr) and the nematic angles Qi{t) are taken in [—7r/2,7r/2). 
For simplicity, we will identify [—vr, tt) with ]R/27rZ and with the unit circle The 
notation di?) describes independent Brownian motions with intensity d cos^(0j —0j) with 
d > 0. We are aiming for a separation of the particles into two (locally dehned) groups, 
one with cos(0i — 0^) > 0 and one with cos(0i — 0^) < 0. Therefore, dihusion must 
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vanish near cos( 0 j — 0 *) = 0 , to avoid particles crossing the boundary between these two 
angular regions. The factor cos^(0i — 0*) is intended to avoid that the two groups mix. 
The hrst term in the velocity evolution describes alignment with the bacteria within the 
interaction radius R> 0 with an alignment frequency z/Sign(cos(0i — 0i)) sin(0i — 0*) 
with u > 0. Due to the factor Sign(cos(0i — ©i)), alignment occurs parallel or anti¬ 
parallel to the mean angle 0j, depending on whether 0j is closer to 0^ or 0j -|- vr. Note 
that Sign(cos(0j — 0*)) sin(0j — 0j) is a vr-periodic function and so, does not depend on 
the choice of 0 j or 0 j -|- tt as the representative of the corresponding angle of lines. 

Remark 2.1 Also other notions of nematic alignment between several particles are pos¬ 
sible. In m the nematic mean direction was defined as 


Qfit) = arg 




Sign(cos(0i(t) - 0fc(t)))n(0fc) 


(2.4) 


The main difference between this definition and (2.3) is that the nematic angle defined 


in (2.3) depends on the particle index i only via the position Xfit), but is independent of 


the direction of the i-th particle, 0*. The outcome of (2.4) on the other hand depends 


on both the position and the direction of the i-th particle. Whilst the definition from B3 
might be closer to reality, if the dispersion around the mean angle is small, both types of 
nematic alignment lead to the same behavior. However, the independence from 0* makes 


(2.3) mathematically much easier to handle. 


2.2 Mean Field Model 

We follow the usual procedure to derive the mean held limit as described for example 
in [23]. The procedure starts with the dehnition of the empirical distribution function 
f^{x, 6, t), with X eMf, 6 eR (modulo 27r) and t > 0 , according to 


where 5(^Xi{t),ei{t)){.x, 6) is the Dirac delta in {x, 6) space located at {Xfif), Ofit)) at time t. 
In the absence of noise {d = 0), f^ is a deterministic measure that satishes a hrst 
order kinetic equation if we assume that all {Xi,Vi) fulhll (2.1)-(2.3). If d 7 ^ 0, f^ is a 
random measure whose law satishes a second order kinetic equation. Then as N ^ 00 , 
f^ approximates a distribution function f{x,6,t), satisfying the following Kolmogorov- 
Fokker-Planck type equation: 


dtf + VoV,, ■ {v{e) f) = Qal/, 
where the collision operator Qai is given by 


Q^]f = i'Sign(cos(^ - 0/)) sm(^ - 0^) / + d cos^(^ - ©/) 5,?/ 


(2.5) 


( 2 . 6 ) 


5 










and where 0/ is given by 


^(20/(a:,t)) (2.7) 

Jf{x,t)=f j v{2e) f{y,e,t)dyA.e. (2.8) 

J\y—x\<R J —TT 


We note from (2.7) that 0/ is an angle of lines, dehned modnlo tt and a representative 
of 0/ will be chosen in [—7r/2,7r/2). We also note that (2.7)-(2.8) can be eqnivalently 
written 

[ f sm2{9 - ef{x,t)) f{y,9,t)dyd9 = 0. 

J \y—x\<R J —TV 


At this point the alignment interaction is non-local which is reflected in the fact that the 
compntations for Jf{x,t) involve integrals over a ball centered aronnd x. As a next step 
we introdnce dimensionless variables. Since we are interested in macroscopic phenomena, 
we nse a hydrodynamic scaling. In the following we will do the non-dimensionalization 
and the hydrodynamic scaling in one step. A more detailed description for a similar 
equation can be fonnd in [21] . On a microscopic scale we nse the reference time to = l/u 
and reference length xq = no to- We also introdnce a scaled diffnsion constant D = dto. 
On the macroscopic scale we want to nse nnits = — and Xq = —, which are large 
compared to the microscopic ones. We introdnce the new (dimensionless, macroscopic) 
variables x = -^ = and t = and set f{x,9,i) = f{x,9,t). 

Additionally we scale R = ^ = = er with r = -^ = 0(1). Note that the scaling of 

R snpposes that the interactions between the particles are local. Another scaling which 
takes better acconnt of the non-locality of the interactions will be investigated later on 
(see Section]^. With this scaling and dropping the hats for better readability we hnd 
that f = satishes the problem; 


e{d,r + V. . (t.(9) /•)) = QIX, 


(2.9) 


with (omitting the variables x and t): 


Qhfi^) = 9e Sign(cos(6' - 0^)) sin(6' - 0^) f{9) + D cos^(6' 


©/) 


def{9) , (2.10) 


and 


v{2e}{x,t)) 


jf{x,t)y 


Jf{x,t) = f j v{29) f{y,9R)dyd9. 

J —TT J\y—x\<er 


By Taylor expansion we get: 


( 2 . 11 ) 

( 2 . 12 ) 


t) — 0j{x, t) + 0{£^), 
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where 


= ( 2 . 13 ) 

jf{x,t)= f v{ 2 e) f{x,e,t)de. ( 2 . 14 ) 

J —TT 

Again, i9j is an angle of lines, defined modulo tt and a representative of 0/ is chosen in 


—7r/2,7r/2). As above, (2.13)-(2.14) can be equivalently written 


/ sin2(6* — 6*j(x, t))/(x, 6*, f) d6* = 0. 

J —TV 

Dropping the terms of order or smaller, we hnd that / = f‘^{x,6,t) satishes the 
perturbation problem: 


Ept/' + V, ■(!.(»)/')) 


QlJ = Do 


Signf cos{0 — Of)') sin(^? — Of) f + D cos'^{0 — Of) def 


where Of is given by (2.13). 


(2.15) 

(2.16) 


3 Properties of the alignment operator Qai 

3.1 Functional setting 

Let K = 1/ D. We can write: 

Sign(cos(0-0/)) sin(0-0/) ^ ^ 

''-cos2(« - Of) -^ 


Q"J = Dde( cm‘(0 - 0,) 


Now, we introduce 


(3.1) 


Definition 3.1 For any angle of lines 9 q G [—7r/2,7r/2), we define the Generalized von 
Mises (GVM) distribution with direction Oq and concentration parameter n as: 


Mooid) = ^exp 


K 


X I cos(6'- ^o)! 

where 

J cos6>0 

is the normalization coefficient such that 


9 e [-IT,it), 


(3.2) 


K 


COS 9 


d9, 


/ Me,i9)d9= / Me,{9)d9 = l. 

Icos9>0 Jcos9<0 


(3.3) 
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Figure 1: Mq{6) for k = 0.5,2,10 (red-dotted, black-solid, blue-dashed). 

Figure shows Mq{9) for different values of k. With this definition, we have 
Lemma 3.2 The alignment operator can he written: 

f 


Qaifi^) = Dde cos {6 - Of) Mg de 


Ml, 


(3.4) 


Proof. Using (3.2) we calculate 

/ 




= dgf - f dg{\ogMg^) 
= dgf + Kf dg 

= dgf + Hf 


cos(0-0/)|. 

Sign(cos(6* — 9f)) sin(6' — 9f) 
cos‘^{9 — 9f) 


which, upon multiplication by cos^(6' — 9f), yields (3.1) and shows the claim. 


We now give a more precise functional setting in which to define the collision operator 
For this purpose, we introduce the following: 

Definition 3.3 (i) Let 9 q G [—7r/2,7r/2) be a given angle of lines. The linear operator 
Cg^^ is defined as 

Cg^if = D dg[cosf{9 - 9o) Mg^ dg(p), (3.5) 

Mg^ 

for any smooth function ip defined on M/(27rZ) with values in M. 

(a) The space Tig^ is defined by 

/ 7T 

|5f(6')|^ Mg^{9) d6' < cxd|. 
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It is a Hilbert space when endowed with the inner product 

i9,h)ng= [ g{9)h{e)Mea{0)de, 


'Hoq ^9 : 9 ')rie(^ ■ 


and associated norm 
(Hi) The space Ve^ is defined by 

= (s'e "Heo, such that \/ MefiO) \ cos{9 - 9o)\d0g{9) e L^(-7r,7r)|, 


where deg is defined in the distributional sense. It is a Hilbert space when endowed with 
the inner product 

i.9. h)ve^ = {g, h)ng^ 

p7T 

+ / (VMgfi9) I cos{9 - 9o) I dgg{9)) (Mgfi9) \ cos{9 - 9o) \ dgh{9)) d9, 


and associated norm Us'llve^ = {9y9)veg- 

(iv) Cg^ defines a continuous mapping from Vg^ into the dual space Vg^ ofVg^ by: 


{^9oT,i^)v'^^,Veo 


= -D 


[\/Mgfi9) I cos(6' - 9o)\dgip) {\/Mgfi9) I cos(6' - 6'o)| dgfi) d9, (3.6) 


for any (p, E Here, (•, •)v^ ,Veg denotes the duality bracket between and its 

dual This duality bracket extends the inner product ofTig^, i.e. for f G Tig^ C Vg^ and 
g G Voq, we have {f,g)v'g ,Veg = / f g^do^^- The operator Cg^^ is formally self adjoint, 
i.e. for all (f, fi E Vg^. 

(v) The operator is defined for functions f E L^(—tt, vr) (which guarantees that 9f is 
defined) such that f /Mg^ E Vg^ by: 

f 


= MgC-g 


f ^f\Mn 


(3.7) 


By abuse of notation, we will simply write (3.6) as follows: 


{Cg^(p, '0)v'^,V9o =-T> dgT de'ip Mgfi9) cos^(0 - 9q) d9. 

J —TV 

Now we investigate the solvability of the equation 


(3.8) 


for a given function fi. In view of (3.6), this problem can be recast into the following 
variational formulation: 


-D 


dg(fi dsx Mg, cos’(9 - Og) d9 = / -i/i x Mg, AO, 'Ox e Vg,. 


(3.9) 

(3.10) 


For this variational formulation, we have the following 
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Proposition 3.4 (i) (Homogeneous problem) Suppose '0 = 0. Then, the unique solutions 


of (3.9), (3.10) are piecewise constant functions, i.e. there exist C± G M such that 


(p{e) = 


C+ if cos{9 — Oq) > 0, 
C_ if cos(0 — Oq) < 0. 


(3,11) 


(a) (Nonhomogeneous problem) Lettp G Then, there exists a solution to (3.9), (3.10) 
if and only if the following two solvability conditions hold simultaneously: 


/ 'ipMgf^d6 = 0, / 'ipM0fjd6 = O. (3-12) 

I cos(9—9q)>0 Jcos{9—9q)<0 


Under these solvability conditions, there exists a unique solution ip to (3.9), (3.10) which 
satisfies the additional conditions 


/ Me^ d9 = 0, ip Mq^ d9 = 0, 

fcos{9—9Q)>0 Jcos{9—9o)<0 


(3.13) 


and the general solution in Ve^ to (3.9), (3.10) is the sum of the unique solution satisfying 
1(3.13) and the general solution (3.11) of the homogeneous problem. 


(in) (Invariance) Let 9 q = 0, 0o ^ Vq satisfying the solvability condition (3.12) and (pQ be 


the unique solution to (3.9), (3.10) satisfying (3.1^ . Let now Oq G [— 7r/2,7r/2), 9o ^ 0 
and define and ip such thatfj{9) = — ip{9) = ipo(9 — 9o). Then, ip is the unique 

solution to (3.9), (3.10) satisfying ( 3.H^ with that value of 9 q. 


Proof. We easily verify that we can choose = 0 and thns we drop the indices Oq in the 
dehnitions of the operators and spaces. 


(i) W e first verify that the fnnction ip given by (3.11) belongs to V and satisfies (3.9), 
(3.10) with '0 = 0. We clearly have ipy/M G L^(—7r,7r) and 


dgip — {C+ — CJ) (5_^/2 — d-K/2)^ 

where 5±,r/2 is the Dirac delta at 6^ = ±7r/2. Since | cos 6^1 a/M is a C°^ fnnction which 
takes the valne 0 at 9 = ±7r/2, we have |cos6*| \fMdgip = 0 G L^(—tt, vr) and ip satis- 
hes (3.10) with right-hand side eqnal to zero. Now, we snppose that 93 G V is a fnnction 
sat isfyin g (3.10) with right-hand side eqnal to zero. Taking x = as a test fnnction 
in (3.10), we get | cos^l \/M deip = 0. This implies that ip is constant on all connected 
components of the set (| cos^l y/M)~^(R.) and is 27r-periodic, which implies that ip is of 


the form (3.11). 


(ii) Let 0 G V be snch that there exists a solntion ip of the variational formnlation 


(3.9), (3.10). Taking y of the form (3.11), we get that 0 mnst satisfy the solvability 


conditions (3.12). Now, we denote by V the snbspace of V dehned by 


V = < (p G V snch that / 

^ J CC 


ip M d9 = 0 and 


ip M d9 = 0 


cos 6>0 


' cos 6<0 
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Suppose that 'ip satisfies (3.12) and that (p is a solution of the following variational 
formulation: 


(pG V, 

D / dgip dex M cos^ 9 d9 = 


'ipxMd9, VxeV. 


(3.14) 

(3.15) 


Then, we claim that 93 is a solution of (3.9), (3.10). Indeed, let x Define X = X~ Xo 


with Xo constructed according to (3.11) with the constants: 


6 *+ = 


xMd9. 


' it cos 6>0 


Then, x Therefore, we can apply (3.15) with test function x- But since Xo D of the 


form (3.11) and since ip satisfies (3.12), the contribution of xo to both the left and right 


hand sides of (3.15) vanish. Consequently, cp satisfies (3.9), (3.10) as claimed. 


We now show that the variational formulation (3.14), ( 3.15[ ) admits a unique solution. For 
this purpose, we apply Lax-Milgram’s Lemma. First, we note that V is a closed subspace 
of V. Indeed, it is obvious that the linear forms ip —)■ M d9 are continuous. 


Therefore, V is a Hilbert space. Then, it is clear that the left-hand side of (3.15) defines 


a symmetric continuous bilinear form a{ip, x) on V and that the right-hand side of (3.15) 
defines a continuous linear form L{x) on V. To apply Lax-Milgram’s Lemma, it remains 
to prove that a is coercive. This amounts to proving a Poincare inequality, i.e. that there 
exists a constant C > 0 such that for all (p G V we have: 


\d,ipf M e Ae >c \ip\'’-Md0. 


(3.16) 


We show that 

[ \deip\‘^ M cos^ 9 d9 > C [ \ip\^Md9. (3.17) 

Jcosd>0 Jcosd>0 

Obviously, if such a result is true, it is also true if the integration domains are changed 


to the set jcos 9 < 0} in both integrals, and this implies (3.16). We now concentrate on 
obtaining (3.17). Let G V and let 6*1 and 02 be such that cos0i > 0, cos 02 > 0. We 
have: 


r-02 


iP{92)-iP{ 9,)= deiP{9)d9. 

Joi 

Multiplying by M(0i), integrating with respect to 0i on (— 7 r/ 2 , 7 r/ 2 ) and using that 
/cos 0 i>o^(^i) ^(^ 1 ) = 0 and get: 

^( 02 ) = [ ( r deiP{9) d 0 ) M( 0 O d 0 i. 

7COS 01 >0 ^^01 ^ 

Taking the square, multiplying by M( 02 ) and integrating with respect to 02 on {—n/2, n/2) 
leads to: 


/ 1 ^( 02 ) pM( 02 )d 02 = / 

/cos 02 >0 -/cos 02 >0 


' deip{9) d0) M(0i) d0i ^ M( 02 ) d 02 . 


'COS 01 >0 ^ ^ 01 


11 
































and by the Cauchy-Schwarz inequality with respect to the inner product on "H, again 
using that /cos 0 i>o ~ obtain 


1^(02) rM(02)d02 < 


' COS 02 >0 


J COS 01 >0 J cos 02 >0 

< Ji + . . . + /4, 


r>6»2 2 

de'tp{e)d6 M(0i)M(02)d0id02 


'01 


(3.18) 


where /i,... /4 are the integrals with respect to the sets Oi to 04 with Oi = {| sin0i| < 
1 /2, I sin6'2| < 1/2}, 02 = {| sin6*i| < 1/2, | sin6^21 > 1/2}, O3 = {| sin^il > 1/2, |sin6'2|< 
1/2}, O 4 = (Isin^il > 1/2, |sin6'2| > 1/2} respectively. For |sin6*| < 1/2, there exist 
constants Ci and C2 such that 0 < Ci < M{ 6 ) < C2 < 00 and 0 < Ci < cos^ 6 M{ 6 ) < 
C2 < 00. So, for (6*1, 6*2) G Oi, using Cauchy-Schwarz inequality, we have: 


r ‘02 


de^PiO) de<C 


'01 


\^e^p{ 9 )\^ M( 9 ) cos^^d^ 


' COS 0 >O 


1/2 


(3.19) 


and consequently, 

h<C [ \^e^l^{9)\‘^M{9) cos^9d9, (3.20) 

J cos 0 >O 

where here and in the following, we denote by C generic constants. Take now {9i, ^ 2 ) G O 4 . 
Then, using Cauchy-Schwarz inequality, we can write: 


-02 


^02 


'01 


'01 


dgij{ 9 )d 9 = I |cos 0 | ^M{ 9 )dei){ 9 ) (|cos 0 | ^/M{ 9 )) d 9 

/ |‘02 J /1 \ 

'01 


< 


^ ' \dei){9)\^ cos^9M{9) d0^ 


cos- 


d 9 Y 
' 9 M{ 9 )) 


Now, we remark that 


and deduce that 
d9 


do 


sin 9 


^02 


cos‘^9M{9) Kjg^ 
1 


do 


M{ 9 ) 

1 ^ 1 


= K 


M{ 9 ) cos^ 9 ’’ 


d 9 


M{ 9 ) J sin 6 * 

1 1 


K, \M(6*2) sin 6*2 M{9i) sin 6*1/ k, Jg_^ M{9) sin^ 6* 


1 

+ - 


I'd 2 


cos 9 


Then, we have 


M{9i)M{92 


-02 


d 9 


Ig^ 9 M{ 9 ) 


1 [M{9i) M{92 


1 

H— 

K 


M{9i)M{92) cos 9 


M{ 9 ) 


Slid 9 


K \ sin 6*2 sin 6*1 
The hrst term is bounded on O 4 . For the second term we remark that 

1 


d 9 . 


M{9i)M{92)/M{9) = exp 


K 


COS 9 


1 1 
cos 6*1 cos 6*2 


(3.21) 


d 9 . (3.22) 


(3.23) 
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Now, 1/ cos^ is a convex function and 9 G [ 6 *i, 6 * 2 ]. Assuming, without loss of generality, 
that 61 < 02 , we have: 


1 02-9 1 9-9i 1 

cos 9 ~ 02 — 9i cos 9i~^ 02 — 9\ cos 02 ’ 


and we deduce that 


1 

cos 9 



9x-9 1 9-02 1 

02 — 9\ cos 9\^ 02 — 0\ cos 02 


< 0 . 


Therefore, M{0i)M{02)/M{0) < 1 and since is also bounded for 0 such that sin 6 ' > 
1 / 2 , we deduce that the second term of (3.23) is also bounded. Consequently, for ( 6 ^ 1 , 02) G 
fl 4 , we get 


(*02 2 

de'ip{9)0i9 M{9i)M{02) <C 


\dei>{0)\^ M{0) cos^0(\.0, 


(3.24) 


'01 


' cosd>0 


which implies, upon integrating with respect to ( 6 ^ 1 , 6 ^ 2 ), that 


l 4 <C [ \^e^p{0)\^ M{0) cos^9d9. (3.25) 

Jcos9>0 

Now, we take ( 6 * 1 , 6 * 2 ) G VL 2 (obviously the case of ( 6 * 1 , 6 * 2 ) G fls can be treated similarly). 
In this case, I 2 is decomposed into I 2 = I 2 + /| where and /| correspond to integrals 
over the domains = (“''’’/S, 7r/6) x (—7r/2, —vr/b) and Vt\ = (—7r/6,7r/6) x {'k/Q,'k/2) 
respectively. We consider J| and a similar proof can be made for J^. For ( 6 * 1 , 6 * 2 ) G 
we can decompose: 


[ ' de^{0) d9= f ^ deij{9) d0 + f ' deij{0) d0, 

^ 0l j 61 j 7T f 6 


and consequently we have: 


-02 


dei!{0) d9 


'01 


< 2 


t/6 


deip{0) d9 


'01 


+ 2 


f*02 2 

dg‘ip{0) d6' . 

t/6 


(3.26) 


Each of the terms will be evaluated separately. For the hrst term, we proceed as for Ji and 
notice that M{0) and M{0) \ cos 6^1 are bounded from above and from below away from 
zero, so that we can estimate it similarly to (3.19) and its contribution to 


/| leads to a 


similar estimate as (3.20). For the second term, we use (3.21) with 0i replaced by tt/Q and 


likewise, we hnd (3.22) still with 6*1 replaced by tt/Q. We can keep on going through all 


the steps of the proof and we eventually get an estimate similar to (3.24) with 0i replaced 


by 7 r/ 6 . This shows that the second term of (3.26) leads to an estimate similar to (3.25). 
Therefore, we get 

l 2 <C [ \^e^p{0)\^ M{0) cos^9d9. (3.27) 


' cos6>0 
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Now, inserting (3.20), (3.25), (3.27) and its counterpart for J 3 into (3.18) leads to the 


Poincare estimate (3.17|) and eventually to (|3.16|), which ends the proof of the Poincare 
estimate. 


Ill 


We restore the indices Oq. Then it is obvious that i/’ and ip satisfy the conditions (3.12), 


the 


(3.13) respectively and that they belong to Ve^- So, they belong to Veg. Furthermore, by 
range of variable 9' = 9 + Oq in the variational formulation satisfied by ipo (^i.e. 


with 


9q = 0), we easily see that (p satisfies the variational formulation (3.9), ( 3.10[ ). Since the 
solution of this variational formulation in the space is unique, ip is this solution. ■ 


3.2 Equilibria 

We now define the notion of equilibria and characterize them: 

Definition 3.5 An equilibrium is a function f: 9 G M/(27rZ) 1 —>■ f{9) such that f G 
L^(—7r,7r), f /Mfj^ G Vq^ and Q^aiif) — 0- denote by S the set of equilibria of 

We also introduce 


Definition 3.6 Let p_ be real numbers in [0, cxd) and 9 G [—7r/2,7r/2) be an angle of 
lines. The function fp^^p_fi is defined by 


f p+Mg{9) for cos{9 — 9) > 0, 
\ P_Mq{ 9) for cos{9 — i9) < 0. 


(3.28) 


It can also be written: 


fP+,P-,9 = Mg{9) [p+ xt + P- Xg ] , 


where ore the indicator functions: 


f 1 if ±cos{9 — 9)>0, 
I 0 if ± cos( 6 '— 0 ) < 0 . 


(3.29) 


(3.30) 


Now, with this definition and Prop. 3^, we can state the 
Lemma 3.7 The set £ of equilibria of is given by 

£ = {fp+,p-,e I (p+,p-) e [0,CX))^ 0G [-7r/2,7r/2)}. 


Proof. We use (3.7). Thanks to Prop. 3.4 (i), we first note that, for a given 9q G 
[— 7 r/ 2 , 7 r/ 2 ), we have 


Meo^9o 


f 


M, 


do 


= 0 


/ . 




9o 


is of the form ( |3.11 ) 


3(p+,p-) G [0,oo)^ s.t. / = 


00 - 
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Applying this statement with Oq = 9f, we deduce that if / G i.e. if Q^ J = 0, then / is 
necessarily written as / = and is consequently of the form (3.28). Conversely, for 

given (p+,p_) G [0, cxd)^, and 6 G [—7r/2,7r/2), we consider fp_^^p_^g and need to show that 
Qaifp+,p-,e = 0- Since fp^^p _^0 is such that /Mq) = 0 (by the equivalence shown 

above), it is enough to show that Of - = 9, or equivalently, that: 

■> p_|_ ,p_ 


/ sin2(0 - 9) /p+,p_,e-d0 = 0. 

J —TT 

But we can make the change of variables 9' = 9 — 9 and use the periodicity to write: 
f sin 2(0 - 0) fp^,p_, 0 d 9 = j sin20/p+,p_,o d0, 

J —TT J —TT 

and this last integral is clearly zero by symmetry. 


3.3 Generalized Collision Invariants 

We hrst dehne the notion of a collision invariant, thanks to the 

Definition 3.8 A function 9 G M/(27rZ) i—)■ d'(0) is a collision invariant (Cl) if and only 
if it satisfies 

/o“.,/(9)«(9)d9 = 0 for all f. 

We deliberately leave the function setting vague as this notion will have to be enlarged 
later on. Clearly, the set of collision invariants is a vector space. At this point, however, 
we only have a one-dimensional space of collision invariants, spanned by the constant 
function d'(0) = 1, corresponding to mass conservation. Following the approach in [23] . 
we therefore generalize the concept of a collision invariant and dehne a generalized collision 
invariant: 

Definition 3.9 Given an angle of lines 9 q G [—7r/2, n/2), a function G is a gener¬ 
alized collision invariant (GCI) associated with 0o if and only if 

(>Ceo^,^)v'o,Veo = 0 for all f e Ve^ A L\-7r, n) s.t. 0/= 0o. (3.31) 

Here the eguality 9f = 9o is an eguality of angles of lines i.e. holds modulo vr. The set of 
GCI associated with 9 q is a linear space denoted by Qe^- 

We note that for / G Ve^ such that Cg^j^ G Heo, the duality bracket {Cg^j^, \k)v' 
reduces to the classical integral / Cggj^{9) dt(9) Mgg{9) d9. Since being a Cl equiva¬ 
lently means that 

/£s,^(9)>l'(9)iVy(9)d9 = 0 forall/, 
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we pass from the definition of a Cl to that of a GCI by replacing the orientation of / by 
an arbitrary orientation 6q and requesting that the integral is zero only for those functions 
/ whose mean orientation 6f coincides with Oq. Because this means we restrict the set of 
functions /, we expect to hnd more 'h, fulhlling the dehnition of a GCI as compared to 
that of a CL We can reformulate the dehnition of a GCI in the following way: 

Lemma 3.10 Given an angle of lines 9o G [—7r/2,7r/2), a function 'L G Vqq is a GCI 
associated with 9o if and only if: 

3a G M such that = asm2{9 — 9o) in (3.32) 


Proof. We note that the constraint 9f = 9^ can be equivalently rephrased as the linear 
constraint: 

/ sm2{9-9o) f{9)d9 = 0. 


By a classical duality argument already used in |23], statement (3.31) is equivalent to the 
existence of a constant a G M such that 

® / sm2{9 - 9o) f{9)d9 for all / G G L^(-7r, tt). 

J—TV 


Using the formal self-adjointness of Cg^ (see Dehnition 3.3 (iv)), this is equivalent to 


= «(sin2(0 - ^o), for all f eVg^G L^(-7r,7r), 


/ 


which implies that d' satishes (3.32). The converse is obvious. 


We now have the following result: 


Proposition 3.11 (i) Given an angle of lines 9o G [— 7r/2,7r/2), the space Qg^ of GCI 
associated with 9 q is of dimension three and is spanned by gg^, with x^^ given 

by (3.30) and gg^, the unigue solution in of the variational formulation (3.9), (3.10) 
with right-hand side fjiO) = sin2(6' — 9q) and satisfying the cancellation condition (3.13). 
(a) The functions Xg^ o,nd gg^ can be written: 


XggiO) = XoiO - 9o), geoiO) = goi9-9o). (3.33) 

In the remainder, we will write x^ 9ind g for Xo ^.nd g^. In addition, g satisfies the 
symmetry conditions: 

g(-0) = -g(e), g{7r-9) =-g{9). (3.34) 


16 














Proof, (i) First, we note that the solution space of Problem (3.32) is the linear space 
spanned by the solution of the same problem with a = 1. Therefore, we restrict ourselves 
to a = 1. Now, solving Problem (3.32) with a = 1 amounts to solving the variational 
formulation (3.9), (3.10) with right-hand side 'ip{0) = sin2(0 — 6o)- We apply Prop. 


3.4 (ii). A necessary and sufficient condition for the existence of a solution is that the 


right-hand side satishes the solvability condition (3.12), which in the present case, with 
■0 = sin 2(0 — 0o) reads: 


' cos(6—(9o)>0 


sin 2(0 — 0o) Moq d0 = 0, 


' cos(6—6o)<0 


sm2{e-9o)Me, d0 = 0. 


This condition is obviously satished by symmetry. Therefore, there exists a unique solution 


qoq of the problem satisfying the cancellation condition (3.13) and all solutions of the 
problems differ 


rom this one by the addition of a solution of the homogeneous problem, 

± 

'o' 

Finally, 


i.e. of the type (3.11). Such solutions are arbitrary linear combinations of Xeg- 
(ii) The statement is obvious for x%- For 9eo: we simply apply Prop. 


3.4 


111 


we remark that sin 20 satishes the symmetry conditions (3.34). Therefore, using the 


uniqueness part of Prop. 3.4 (ii), we deduce that g satishes the same symmetry conditions. 


In the present simple two-dimensional case, we can have an explicit formula for g, 
given by the Lemma below and depicted in Figure 


Lemma 3.12 For 9 G [0,7r/2] we have 

TT 

sin2Q; exp(- 


9(0) = - 


-) da 


cos2/3 exp(-^) 


d/9, 


and g is extended to ]R/(27rZ) by the symmetry conditions (3.34)- 

Proof. Solving for g consists of solving the following problem: 

d / 9 k; do / \ ^ 

— (cos 0e — (0)) = sm20 e 
d0 V d0 / 


(3.35) 


(3.36) 


Eq. (3.36) can be solved explicitly for 0 G [0,7r/2] and then extended using the sym¬ 


tions (3.13). 


metry conditions (3.34). Integration constants are determined by the cancellation condi¬ 


4 Hydrodynamic Model 


4.1 Derivation 


This section is concerned with the limit £ —)■ 0 in Problem (2.15)-(2.16). The hrst step 
is the observation that the dominant term Q^if^ = 0 requires the limiting distribution to 
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- 1-1 


Figure 2: g{9) for k = 0.5,2,10 (red-dotted, black-solid, blue-dashed). 


lie in the kernel of Qau characterized by Lemma 3.7 To determine the equations satished 


by the macroscopic quantities p+(x,t),p_(a:, f) and 6{x,t), it is necessary to consider the 


terms of order £ in (2.15). This is done by multiplication with the three independent GCI 


given in Prop. |3.11 
those in 


and subsequent integration over 9. The computations are similar to 


and are described in more detail in the Appendix (Section 8.1). Here, we 
simply state the result: 


Theorem 4.1 Taking the (formal) limits —)■ 0 in (2.15)-(2.16), we obtain 

where the macroscopic quantities p±{x,t) and9{x,t) have values in [0, cx)) and [—7r/2,7r/2) 
respectively and fulfill the following system of equations: 

dtp+ + diW^- {p+ v(9)) = 0, (4.37) 

dtp- - di Va: • {p-v{9)) = 0, (4.38) 

{p+ + P-) dt9 + d 2 ip+ - P-) {v{9) ■ Va,)9 + pv{9)-^ ■ Va:(p+ - p-) = 0, (4.39) 

where we recall that v{9) = (cos 0, sin 0) and v^{9) = {—sin 9, cos 9). The coefficients di, 

d 2 and p are given by: 


di = (cos) 


M, 


do — 


{ 9 ^ s)m ^ 1 (f/sin)M 

(0>m’ 


(4.40) 


where g is the GCI defined in Prop. 3.11 and explicitly given by (3.35). For a function 
(p{9), we denote by {(p)m the following average: 


^7r/2 2 

— 2 / ip{9) AI{0) do — — / do. 


(4.41) 


■‘K JQ 


where M is My for 9 = 0. 
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The proof can be found in the Appendix (Section |8.1[ ). 

As shown in the proof of Lemma 3.7, the macroscopic quantity 6{x, t) is the local mean 
angle of lines as dehned in (2.13) of the equilibrium distribution /. The functions p+{x, t) 
and p_ {x, t) describe the local densities of particles that point in dir ectio n \0{x ,t) — 
Tr/2, 9{x,t) + ti/2) and \9{x,t) + 7i/2,9{x,t) + 37r/2) respectively. Eqs. (4.37) and (4.38) 
simply describe number density transport in the direction of v{9) (for p_|_) or —v{9) (for 
p_). These equations express local conservations of the number densities of the particles 
moving in the direction of v{9) or —v{9) respectively. Since GCI do not express local 
conservations (unless they are Cl), in general one cannot expect that they will lead to 
conservation laws. Here, Eqs. (4.37) and (4.38) are conservation laws in spite of the fact 
that they are derived from using the GCI This unusual circumstance is due to the 
particular expression of the collision operator 

Note that whenever = p_, it follows that dt9 = 0 and dt{p+ + P-) = 0. An 
alternative formulation of System (4.37)-(4.39) can be obtained by dehning p = + p_ 

and 5 = p+ — p_. We get: 


dtp + d,V^-{5v{9))=f), 

(4.42) 

o' 

l> 

1 

(4.43) 

pdt9 + d 2 6 {v{9) ■V^)9 + p {v{9)^ ■ V.)h = 0. 

(4.44) 


It is instructive to compare with the Self-Organized Hydrodynamic (SOH) system 
derived in [23j, which, in dimension two, takes the form: 

dtp + diV.-{pv{9))=0, (4.45) 

p (^dt9 + d 2 {v{9) ■ V.)0) + p{v{9)^ ■ V.)p = 0, (4.46) 


where the constants di,d 2 and p depend on k and the collision invariants of the corre¬ 
sponding alignment operator. This model is based on polar alignment. Here by contrast 
to System (4.37)-(4.39), 9 is an angle of vectors, i.e. 9 G M/(27rZ) and the particles are 
distributed according to a classical von Mises-Fisher distribution (i.e. their angular dis¬ 
tribution is proportional to and they move (with a certain dispersion described 

by k) in the direction of 9. 


The fact the nematic alignment leads to an equilibrium angle of lines as opposed to 
an angle of vectors is reflected in the existence of two opposing local maxima at 9 and 
9 + 71 of the GVM Mg{9), whilst the classical von Mises-Fisher distribution has only one. 
The factor cos^{9 — 9f) in front of the diffusion operator in Q^i allows for this separation 
into two groups and causes Mq{9) = 0 for 0 = 0 ± 7r/2, i.e. in equilibrium there are 
no particles moving perpendicular to the equilibrium angle of lines. On the other hand, 
both the von Mises-Fisher distribution and the GVM share their general ” bump-like” 
shape and the convergence to a Dirac delta for k —)• oo. At the level of the macroscopic 
equations further similarities and differences become evident. In the nematic case both 
groups of particles have their own conservation law, consequently we have three instead of 
two macroscopic equations. The general structure remains the same: Transport equations 
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for the density complemented by an equation describing the angle evolution, consisting of 
a time derivative, a convection term and a pressure term. However for the nematic SOH 
model both the convection and the pressure term depend on the local density difference 
6. In the case where there is only one group of particles, e.g. p+ = p and p_ = 0 system 
(4.37)-(4.39) reduces to (4.45)-( 4.46[ ) (with different constants), since then the nematic 
character of the collisions does not play a role. 


4.2 Hyperbolicity of the Macroscopic Model 


An important property of the original SOH system (4.45 )-(4.46) is its hyperbolicity. In 


a variant of the SOH model was examined. It was found that in some regions hyperbolicity 
is lost, which affects the behavior of the solutions. Therefore, we examine the hyperbolicity 


of the macroscopic model (4.37)-(4.39). To simplify the analysis we rescale the system 


with t —)■ dit. In the following we therefore abbreviate ^2 = ^ 2/^1 and p = p/di 


Lemma 4.2 The macroscopic system (4-3'7)-(4-39) is hyperbolic 


Proof. We use the alternative formulation (4.42)-(4.44), which leads to simpler formulas. 


we dehne u{x,t) = {p{x,t),6{x,t),9{x,t)) with x = {xi,X 2 ) G yielding 

dtu + Adx^u + Bdx2U = 0 , 


where 


0 

A{p,6,e) = ( cos^ 
0 


cos 6 —6 sin6'' 

0 —p sin 9 

sin 9 ^ cos 9, 


B{p,5,9) = A{p,5,9 -ti/2). 


To determine the hyperbolicity of the system, it suffices to show that the matrix A has 
real and distinct eigenvalues for all 9 G [—7r/2,7r/2). The characteristic polynomial of A, 
charPoly^( 2 ;) is of third degree and has the shape 


s s 

charPoly^(2;) = - cos 9 d 2 Z^ — [jl sin^ 9 + cos^ 9) z - cos 9 (^p sin^ 9 — d 2 cos^ 9^ . 


To understand its roots, we use the following definition of a discriminant of a third order 
polynomial: For p{z) = a^z^ + 02 ^^ + aiz + Oq we define 

discrp = — 27 alal + 1803000102 + O2O1 — do^Oo — 401O3. 

Next we observe that the discriminant of the characteristic polynomial of A, which we 
define as A = discrcharPoiy^, can be interpreted as a quadratic polynomial in X = {5/p)"^ 
with coefficients dependent on c = cos^ 9, namely: 


A(X) = a{c)X^ + /3(c)X + 7 (c), 
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where 


a{c) = 4:dlc^ c{d 2 + fi)—fi , 

/3{c) = c (dlff — 18 d 2 fi^ — 27/i^ — 36 d 2 fi — 8 (I 2 — 20 d\pi\ 


+ 2/i y—dljl + 18d2 + 10 dl + 27jl + 18/1 d 2 j c — /i^(18 d 2 — d\ + 27) 

7 (c) = 4[c(l -/i) +/i]^ 

If A(X) > 0 for all c G [0,1] and all X G [0,1] we have three distinct, real eigenvalues 
and the system is (strictly) hyperbolic. The roots of A(X) are given by 


^1,2(C) = 


-(d{c) ± v'/3(c)2 -4a(c)7(c) 


2 a(c) 


Simple calculations show that for c > 


9(/l+d2) 
9/i+/i 6^2+86^2 


the discriminant A(X) has no real roots 


and is positive for all X G [0,1]. Otherwise we can distinguish between c > s := and 
0 < c < s. In the first case both roots Xi^ 2 (c) are negative are therefore not relevant. In 
the second case we lose hyperbolicity for X > X 2 {c). However a simple analysis shows 
that X 2 (c) > 1 for all 0 < c < s. In fact X 2 {c) —>■ 00 for c —)■ 0+ and c —>■ s“ and 
mincg[o,s] X 2 {c) = 1. Since X G [0, 1] this shows the hyperbolicity of the system. ■ 


5 Weakly non-local scaling 

In this section we want to investigate how the model changes, if we allow for a relatively 
larger interaction radius. With the scaling choice below, this leads to interactions that are 
no longer local in the limit, but rather weakly non-local. This and other scaling choices 
for different models o f sel f-propelled particles have been examined in [2T]. In the scaling 

we used R = 0{e). If we use R = 0{y/e) instead, namely 
;he following expression for J|: 


2.2 


performed in Section 
R = y/er, we are led to t 


Jf{x,t)= / v{2e) f{y,e,t)dyde, 

J —IT J \y—x\<y/e r 


(5.47) 


and ©J is still given by (2.11). Then, we need to consider higher order terms in the Taylor 
expansion of Jj and 0^. I’hese are given by the following: 


Lemma 5.1 We have 


Jf{x,t) = enr^ {jf{x,t)+ej}{x,t) + 

1 r 


jj{xR) = kAJf{xR), k = 


dvrr^ 


'\i\<r 


iei^de = -. 


0 - =ef + eO) + o{e^), e] = • v^{2ef)), 


(5.48) 

(5.49) 

(5.50) 


and where jf, 9f are given by (2.14) (2.13) respectively 
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Proof. We change variables to y = x + y/e^ and get 


Jf{x,t)=e / v{2e) f{x + ^/e^,e,t)d^de. 

J J —TV 

By Taylor expansion of / with respect to ^/e and using cancellations of the order ^/e 


terms due to symmetry, we hnd (5.48), (5.49). Therefore, we get 


= v{29f) + + 0{e‘^). 


Expanding 0j according to the hrst equation of (5.50), we get 

Jf 


\J 


f\ 


= v(2Qf) = v(26f) + 2e6jv-^(26f) + 0(e‘^). 


(5.51) 


(5.52) 


(5.50). 


Identifying the coefficients of £ in (5.51) and in (5.52), we get the second equation of 


From the expansion given in Lemma 5.1, we deduce an expansion of the operators 


Qiif- recall that their expressions are given by (2.10). This expansion is given in the 

Lemma 5.2 We have 


Qlif = Qlif + eQlif + 0{e^), 


(5.53) 


where is given by (2.16) and Q\J is given by 

Q\if = 0] de 


2(5(cos(6' — Of)) sin^(6* — df) — \ cos(6* — 6f)\'^f 


+2D cos{9 — 9f) sin(6* — 9f) dgf 


(5.54) 


where (5(cos(6* — 0/)) is the Dirac delta at cos(6* — 9f) = 0. 


Proof. Formula (5.53) follows from the hrst equation of (5.50) as well as 


e - ^ = 9 - 9f - e9] + Oie"^). 


(5.54). Formula (5.53) is obvious 


Then, inserting this expansion in (2.10) and using elementary calculus, we are led to 


From the previous two lemmas, in the new scaling, the mean-held kinetic equation (2.9) 
is written: 

e(ai/' + V, ■ {v(0) r) - Ql,r) = Qlr + 0 (£"), (5.55) 

and we will now drop the 0(e'^) term. The main theorem of this section is the following. 
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Theorem 5.3 Taking the (formal) limit £ —)■ 0 m ( 5.55[ ), we obtain 

f t fp_^_(^x,t),p-{x,t),0(x,t)i.^') y 

with fp^^p_fi given by {3.28). The macroscopic guantities p± and 9 fulfill 

dtp+ + diVx- {p+ v( 9 )) = 0 , 
dtP- -diVx- (p-v( 9 )) = 0 , 

(p++ P-) dt 9 + d2 (p+ - P-) {v{ 9 ) ■ Vx )9 + pv{ 9 )-^ ■ Vxip+ - P-) 

= 2k i^2 -\- d^ — p — —^ Vx ■ ((p+ + p-)'Vx9 ), 


where di, d 2 and p are given by (f.fO) and 


4 = 2 


“ cos^ / M 


The definition of {.)m can be found in (4.41). 


(5.56) 

(5.57) 

(5.58) 

(5.59) 


The proof of this theorem is given in the Appendix (see Section 8.3). The weak non¬ 
locality of the interactions leads to the addition on a diffusion term on the right-hand side 
compared to purely local interactions. The sign of the diffusion constant depends on the 
sign of V = d 2 + ds — p — Lemma 5.4 shows that P is positive. 

Lemma 5.4 Let the diffusion constant T> be defined by 


T) d2 -\- — p -, 

K 


where 4 o,nd p are given by (4.40) and 4 given by (5.59). Then it holds that 


V >0 for K >0. 


Proof. We start by rewriting 

27 = 4 + 4 — h 


2 

K 


+ 2 


(9^)n.r ^( 9 ^) 

cos^ / M cos^ / 


_(«Mm I o 1(9Sin) 

(9S>m 

1 

WS) 


M 


M 


sm sm 
K -h 2k - 5 


M 


sin —2- 


cos 


COS'" 


2 

K 

sin 

cos^ 


M 


A 

B' 


(5.60) 


Next, using the definition of (.)m given in (4.41), we examine the numerator A. Note that 
since would appear in both A and B, we can omit it. 


A=(9 


sm 

cos^ 


[k (cos^ -1-2 sin^) — cos (cos^ -|-2)] \ 


/ M 


/•^/2 a 

/ [k (cos^ 0-|-2sin^ 6*) — cos 6* (cos^ 6*-|- 2 )] d6* 

I Q COS*^ u 

[-.(2 -y)-,(/ + 2)] d,. 

'0 9 


(5.61) 
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In the last line we used the change of variable y = cos 6 and the dehnition 


G{y) = 5 '(arccos(|/)). 


The last line (5.61) can be split into 
^^G{y) d 


A = 


Jo y ^y 

= Ai + A2. 


(e (2 - y^) dy - 


y {y^ + 2) dy 


y 


(5.62) 


We start by rewriting Ai using integration by parts. This yields 


A, = ^e-‘. {2 - y^) 


y 




2-r 


_^2-|/2 d 
e y -— 

) y ^y 


G{y)dy- / e-yG{y) — 


dy 


d [2-y^ 


<^y \ y 


dy 




where the second equality is obtained by noting that the boundary term vanishes. Note 
that since < 0 on [0,7r/2], which can be seen from its dehnition in (3.35), it also holds 
that 


-^G{y) = -^=L=-^(arccos(?/)) >0 for ye [0,1]. 


(5.63) 


We continue by expanding A\\ 


A\ = / e-yG{y) 




which shows that in fact A\ = —A 2 as dehned in (5.62). Using this and (5.63) we can 
conclude 


A = A} = 


e I- — y--^G{y)dy < 0. 
) y dy 


Finally, since < 0 on [0,7r/2], we observe that the denominator B as dehned in (5.60) 
(again omitting Z^) fulhlls 


5 = = 1^1 g{e)e' 

\ COS^ / M Jq 

which shows the claim and hnishes the proof. 


^ sin^ 

cose- ({0 < 0, 


cos^ 9 
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6 Example: Myxobacteria 

In this section we want to give an example of how to use and adapt the derived nematic 
SOH model to a particular biological problem, the movement of myxobacteria. This 
type of rod-shaped bacteria can glide without the aid of flagella on surfaces. Under 
starvation conditions, prior to the formation of fruiting bodies, they go through a so- 
called rippling phase. During this stage traveling waves of high density bands are observed 
that seemingly pass through each other unaffected. However upon close examination it 
has been found that a large number of bacteria at the interface of two colliding waves 
in fact reverse their direction. These reversals can happen spontaneously, but a higher 
number of head-to-head contacts with other, oppositely moving myxobacteria increases 
their likelihood [321 03] • Biologically these reversals are quite fascinating, because they 
happen by the reorganization of the internal movement machinery as opposed to an actual 
turning process. The interaction of two bacteria that don’t reverse can be described 
by nematic alignment. In terms of other communication between bacteria no diffusible 
chemotactic signal is known to coordinate their movement. At this point the nematic 
SOH model can describe the (nematic) alignment of two colliding bacteria, however the 
reversals still need to be included. 


6.1 Particle and Mean Field Model 

Particle Model 

We complete the description of the particle model from Subsection |2.1| by assuming that 
particles reverse their direction of motion with a frequency, dependent on the density of 
oppositely moving bacteria So we let 0* be changed into 0* -|- tt at Poisson random 
times with frequency Kp-signcos{e^-ep)^ 

p\ := Card{/c G J\f I \Xk — XA < R and ± cos(0fc — 0^) > 0}, 
where Card stands for the cardinal of a set. In other words, we have 

Prob{0i(f + T) = Qi{t) + 71 I 0j(f -|- r) = 0j(f) Vr with 0 < r < T} 

= 1 - exp ( - X{plsigncos(e,-ep) ^) • (6-64) 

(assuming that P*_gigncos (0 -e ) ~ P-Signcos (0 constant in the time interval 

[t,t -|- T]). In these formula, the function A(p±) is a reversal frequency. The way it de¬ 
pends on p± will be discussed later. 


Mean Field Model 

The above considerations change equation (2.5) to 

dtf + VqVx ■ {v{6) f) = Qsxf + Qrevf, 


where the collision operators Qai is given by (|2.6|) and 
Qrevf ix,9,t) = 


fA(cr/,+) f{x,e + ttR) - A(a/_) f{x,6R), for cos(6' - 0/) > 0, 
9 + 7r,t)- A(a/,+) f{x, 9, t), for cos{9 - Qf) < 0, 


(6.65) 


( 6 . 66 ) 
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and where 0/ is given by (2.7)-(2.8) and af^± given by 

1 r r 




RH 


di cos(0—©j^)>0 \x—y\<R 


f{y,G,t) dy de. 


(6.67) 


To obtain local interactions, we assnme as in Section 2^ an interaction radins fnlfilling R = 
0(e) and choose for the reversal freqnency the scaling = toA((j/,±) = 

The local density is then given by 






' it cos(6—0p>O J\x—y\<er 


fiy,^^) dy d0, 


and it can be shown that 


= pf,±{xR) + 0{e‘^), 


where 


Pf,±{^^t) = / _ f{x,o,t)de. 

J di cos{6—0f)>0 


This changes eqnation (2.15) into 


e{a,r + V, . (v{0) /')) = Qhr + eQ“„r, 

where Q^i/ given by (|2.16 ) and 


( 6 . 68 ) 


(6.69) 


Qve^f{x,9,t) = 


KPf,+) fix, 9 + 7i,t)- A(p/,_) f{x, 9, t) for cos(0 - 9f) > 0, 
HPf,-) fix, 9 + 7r,t)- A(p/,+) f{x, 9, t) for cos(0 - 9f) < 0, 


(6.70) 


where 9f and pf^± are given by (2.13) and (6.68) respectively. 


Remark 6.1 In (6.69) the contributions from appear at the same scale as the trans¬ 


port operator, which is a consequence of the scaling choice for the reversal rate A. Using 
Xio'f,±) = toA(cr/,±) instead, would lead to a different model, in which and would 
he of the same order. Consequently also the shape of the equilibria and the macroscopic 
model would look very different. We leave this path for future work. 


6.2 Hydrodynamic Limits 

Macroscopic Model - Local Scaling 


Proposition 6.1 Taking the (formal) limit e: — t 0 m (6.69)-(6.70), we obtain 

f iXydyt'j y fp^(^x,t),p—{x,t),9{x,t)i^'). 
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where the macroscopic quantities p±{x,t) and9{x,t) have values in [0, cx)) and [—7r/2,7r/2) 
respectively and fulfill the following system of equations: 


dtp+ + rfi ■ (p+ v{9)) = A(p+) p_ - A(p_) P+, (6.71) 

dtp- -diV^- {p-v{9)) = A(p_) p+ - X{p+) P-, (6.72) 

{p+ + P-) dt9 + d 2 {p+ - P-) {v{9) ■ Vfi)9 + pv{9)^ ■ Vx{p+ - P-) = 0, (6.73) 

where v{9) = (cos 0, sin 0) and v^{9) = (—sin0, cosi9). The coefficients di, d 2 and p are 

defined in (4.40). 


The proof consists of some simple additions to the proof of Theorem 4T and is given 
in the Appendix (see Section 8.2). 


Macroscopic Model - Weakly non-local Scaling 

To obtain the hydrodynamic limit of the myxobacteria model for the weakly non-local 
scaling described in Section we note that in this scaling 

al+{x,t) = ^^[ _ [ f{y,9,t)dyd9. 

er J±cos(9—0^j:)>O J\x—y\<+:£r 


It holds that 


a 


f^+{x,t) = pf^+{x,t) + 0{e). 


We observe that (5.55) is changed into 


fidj’ + V, . {v(e) /•) - Qi,f - Qi,) = (?»,/■ + 0(e\ (6.74) 


with QIj and given by (5.54) and (6.70) respectively. The resnlting limit is snmma- 
rized in 


Proposition 6.2 Taking the (formal) limit £ —>■ 0 m (6.69)-(6.70), we obtain 

f (^) ^) f) t f1 

With fp^+_p given by The macroscopic quantities p± and 9 fulfill 

dtp+ + di Va; • (p+ v{9)) = X{p+) p- - A(p_) P+, 

dtp- - di • {p-v{9)) = X{p-) p+ - X{p+) P-, 

{p++ P-) dt 9 + d 2 {p+ - P-) {v{ 9 ) ■ Vfij 9 + pv{ 9 )^ ■ Vx(p+ - p-) 

= 2k (^2 -\- d^ — p — —^ Vx • ((p+ + p-)'Vx9), 


can be found in (f-fl) 


where di, d 2 and p are given by (f-fO) and d^, is given by (5.59). The definition of {.)m 
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The proof is a simple combination of the proofs of Theorem 5.3 and Proposition 6.1 


One can observe that the additional reversal term leads to a reaction term on the 


right-hand side of equations (6.71) and (6.72), with the precise shape depending on the 


choice of A(p). Motivated by biological hndings we expect A(0) = Aq > 0 to account for 
spontaneous reversals in absence of other bacteria. Since head-to-head contacts increase 
the reversal frequency, A(p) should be an increasing function of p. To demonstrate what 
type of dynamics could be expected, we give an example, where we choose 


A(p) — Aip^ -|- Aq, 


(6.75) 


with Aq, Ai > 0. For a constant direction, 6{x, t) = 9 and moving along the characteristics, 
the corresponding phase portrait is depicted in Figure]^ One can observe that wherever 
the total density -|- p_ > one of the two groups dominates the other one, whilst 

for smaller values of -|- p_ both attain the same number. This is similar to the waves 
observed in myxobacteria, where at the high density peaks of the waves, almost all bacteria 
move in the same directions, whilst in between the waves, there can be small numbers of 
bacteria moving in different directions. 
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Figure 3: Local dynamics for X{p) given by (6.75). The arrows mark the flow held in 
the (p+, p_)-plane. The red-dotted and green-dashed lines show the values for which 

A(p+)p_ — A(p_)p+ = 0. The blue-solid line shows the threshold values p^ + P- = 2w 


A detailed analytical and numerical analysis for this and other choices of X{p) will be 
examined in a forthcoming work. 









7 Conclusion 


In this work we have derived macroscopic equations for self-propelled, rod shaped parti¬ 
cles moving in two space dimensions starting at the particle level. The type of interaction 
considered is nematic alignment of polar particles, which is commonly found in biological 
and physical systems, but systematically derived hydrodynamic equations for this type 
of alignment are largely missing from the literature. In order to take the limit in the hy¬ 
drodynamic scaling, we used the concept of generalized collision invariants. Two types of 
scalings have been examined: a local scaling and a weakly non-local scaling. The result¬ 
ing macroscopic equations share some similarities with their polar alignment counterparts, 
but differ signihcantly in other aspects. A key difference is the splitting of the particles 
into two groups moving into opposite directions. Lastly, including density dependent re¬ 
versals, we have demonstrated how the equations can be used to describe the behavior of 
myxobacteria. The research in the area of pattern formation and non-equilibrium states 
of large assembles of particles is very active and more macroscopic models are needed to 
better understand the phase transitions and underlying causes. Several extensions of this 
work are possible: to describe a wider range of phenomena it could be generalized to three 
and more space dimensions. Further at this point the limits examined are on a formal level 
and lack a rigorous treatment. Extensive numerical experiments need to be conducted to 
compare the particle model with the macroscopic models and to characterize the different 
phases and states. To better understand the rippling behavior of myxobacteria different 
choices of reversal rates can be examined, both numerically and analytically. 


8 Appendix: proofs 

8.1 Proof of Theorem 14.11 

8.1.1 Formal limit £ —)■ 0 


All limits taken are formal. We suppose that converges to some /° when £ —)• 0. From 
) we see that the limiting distribution has to fulhll Qai/° = 0 and therefore, use of 
Lemma 13.71 leads us to 





(8.76) 


To derive the equations determining p+(x,t), 9{ x, t), we multiply (2.15|) by the 


three GCI xf S'lid associated to given by Prop, 
to 9. '' 


3.11 


and integrate with respect 


8.1.2 Use of GCI xt 

We interpret the integral of against xf^^ follows: 


Ql,n9)xtMd9={C5^ 




Mn 




m 






29 











and, by Definition 3.9, this integral equals zero. Upon division by e we are left with 


' dl cos{0 — 9fe)>0 


{dtr + v{e)-vj^) d0 = o. 


Letting e —)■ 0, nsing (8.76) and the fact that 9fe —)■ 0, we are led to 


{pt + v{9) ■ Vx) fp+{x,t),p.{x,t),e(x,t){(^) — 0. 


(8.77) 


J it cos(0—6)>0 

Note that the integration domain depends on x and t. We calculate 


{dt + v{e) ■ Vx)fp^,p_pe) = M-pe) {dt + v{e) ■ Vx)p: 






sin(6' — 6) {pt + v{6) ■ Vx)d for ± cos(6* — 0) > 0, 


cos^(6* — 6) 
which can also be expressed as: 

{dt + v{e) ■ = Me{9) [xl{dt + v{9) ■ V.)p+ 

+Xs {dt + v{9) ■ Vx)p- + ^ {xbP+ - XjP-) {dt + v{9) -Vx)^ . (8.78) 


We note that in spite of the presence of indicator functions in (3.29), is continuous 

(and even C°^) across the points 9 where cos(6' — 0) =0, so that there is no delta 
distribution appearing at these points when we differentiate it. In the following we perform 
the calculations for the integral over 9 such that cos(6' — 0) > 0, the case of cos(0 — 0) < 0 


being treated in the same manner. The left-hand side of (8.77), denoted by L, is then 
equal to 


L = 


' cos(6—6)>0 


M-P9) d9)dtp+ 


Mg{9)v{9)m-V^p+ 


' cos(6—6)>0 
+ Kp+{ 


KP+ 


/ ) d,e 

lcos{e-e)>o cosp9-9) 

[ Ms{9) 


/cos(0-d)>O cos^(9-9) 

= Li L 4 . 


(8.79) 


Thanks to the normalization condition (3.3), we have Li = dtp+. Using that 


v{9) = cos{9 — 9) v(9) -h sin(6' — 9) v~‘~(9), 


(8.80) 
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together with cancellations due to symmetries, L 2 simplihes into 


L 2 = / M 0 ( 6 *) (cos( 6 '— 6 >) n( 6 *) ■ Vx + sin( 6 * — ■ Va;)p+ 

J cos(6—6)>0 


/ _ Mq{0) cos{6 - 9) d6'j {y{9) ■ Vx)p+ 

'cos(0—6)>0 ' 

M(9) cos6'd6'') (^(19) ■ Vx)p+ 


cos6>0 

= di {v{e) ■ Va:)p+, (8.81) 

where di is dehned in ( 4.40[ ). We easily get L 3 = 0 by symmetry. Finally, using again 
(8.80), cancellations due to symmetry and integration by parts, we hnd that 


^4 =^P+[ M-,{e) 

/cos(0-0)>o cos^e-e) 


= -p+ 

= p+ 


[ _ sm{e - e)de') {v^{e) ■Vx)9 

' cos{9—d)>0 

M-q[9) cos(0 - 9) d^) (n^(0) ■ 


' cos(6—0)>O 

= di p+ (n-^(i9) ■ Va;)d. 


(8.82) 


Finally, collecting the contributions from above, we get (4.37) and (4.38) 


8.1.3 Use of GCI 

We start as in the previous section. We interpret the integral of Q\\f^ against g^^^ as 

/ (?;,/'(«) 3i,. W) id = Ui,. m,. (»)),„ „ . 

and again by Dehnition 3.9 this integral vanishes. Upon division by £ we get 

j {a,r + 1.(9) ■ V.f) 95 ,, d9 = 0. 

Again, we assume that —)■ fp_^_(x,t),p-{x,t), 0 (x,t) 9fe —)■ 9 and hnd 


0 = 


ged9 = X. 


(9, + t.(9)-V44, 

We use the calculations in (8.78) and split X into X = Xi + ... + X 3 with 
^1 = 


{xt 9tp+ + Xe dtp-) ggd9, 


X 2 — 


^(^) ■ {xt VxP+ + Xg ^xP-) + 


sm{9-9) , ^ - \ao\ AO 

K— P+-Xg P-)dt9 ggd9, 


X3 — K, 


cos^{9 — 9) 

cos2(^g -1) d+ - Xg- P-) W) ■ Vx)9 g-g d9. 
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For the first integral inside the expression of Xi we have 

I xt9ed0 = 0, 


using the symmetry relations (3.34) and therefore Xi = 0. For X 2 , using (8.80) and 


cancellations due to symmetry, we get: 


X, = / 


6{x,t) 


sm 


{6 - 6) v^{e) ■ {xg V^P+ + Xg ^xp-) + 
sin( 6 * — 9) 


K 


cos^{6 — 9) 


{Xs P+- Xg P-) dt9 ggd9. 


Now, it turns out that the integrals with respect to 9 on the two sets cos( 6 ' — 9) > 0 and 
cos( 6 ' — (9) < 0 are opposite by symmetry and we get 


X 2 = (gsin}MV-^(9) ■ V^(p+ - p-) + k (g 


sm 


COS ^/M 


{p+ + P-)dt9. 


Finally for X 3 we again use (|8.80|) and cancellations due to symmetry and get: 
X 3 = K Mg(, 


S(x,t) ^- 7 ? — S {xt P+ - Xg P-) {v{0) ■ ^x)9 gg d0 


cos( 6 * — 9) 

= '‘(3—) (.p+-p-){v(S)-v,)s, 

\ cos/m ^ ^ 

since now, the integrals with respect to 9 on the two sets cos{9 — 9) > 0 and cos{9 — 9) < 0 
are equal by symmetry. Combining these equalities, we get 

0 = «(^) + {p^-p.){vi9)-X.)9 

\ COS"' / M \ cos / M ^ ' 

+ (5(sin)MU-^(0) ■ V^(p+- P-). (8.83) 


Dividing this equation by k we get (4.39) with the coefficients dehned by (4.40). 


8.2 Proof of Proposition 6.1 


Proof. The proof consists of some straightforward extensions to the proof of Theorem 
4.1 First of all we note that the limiting distribution still has to fulhll Qai/° ~ 0- 


Therefore we get that 


f fp^{x,t),p—{x,t),d{x,t)i.^^ 


and we are left with determining the integral of (6.70) against the generalized collision 
invariants xf ^ and gg^^. We start by noting that 


QLfp^,p.,m = 


(A(p+) p_ - A(p_) p+) Mg{9) for cos( 6 ' - 6 >) > 0, 
(A(p_) — A(p+) P-) Mg{9) for cos( 6 ' — i9) < 0. 


(8.84) 
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With this it follows easily that 


Qrevfp+{x,t),p-(x,t),d(x,t)Xg^^^ i('^(P+)P— -^(P—)P+)- 


(8.85) 


Finally following the same arguments as in Section 8.1.3 about Xi, we hnd that 


Qrev fp^{x,t),p—(x,t),9{x,t) 99 0) 


which hnishes the proof. 


8.3 Weakly non-local scaling: proof of Theorem |5.3 

Compared to the proof of Theorem |4.1[ we need to compute the contributions of 


( 8 . 86 ) 


We hrst compute: 




P+.P-> 


= P+ 


' cos{9—9)>0 


v{29)Mg{9) d9 + p_ 


' cos(9—9)<0 


v{29)M-g{9) dP. 


Thanks to 


v{29) = cos {2{9 — 9)) v(29) + sin (2(6* — 6*)) ^■‘■(26*), 
and owing to cancellations by symmetry, we get 


= (P+ + P-) { I cos(29)M(9) d9 j v(29) 


' cos 6>0 

= ip+ +P-){cos{2-))mv{29), 


and consequently 

b/p+.p_,el = (P+ + P-) (cOs(2-))m, 

since (cos(2-))m > 0. Now we get: 

= k {cos{2-))m + p_)v{29)) 

= k {cos{2-))m (^^xip+ + P-) v{29) + 4Va;(p+ + p-) ■ Xj v-^{29) 

+4(p+ + p_) {Ajv^{29) -\Vj\^v{29))y 
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and finally 


1 


6} 

J p+,p-,e 


(d. 




2I7V -I ^ •'p+,P_, 0 ^ 

l-'/p+,P_,el 


2k 


P+ + P- 
2k 

P+ + P- 


(Va;(p+ + P-) ■ 'VxG + (p+ + P-) 
X ■ ((p+ + P-)Vx^)- 


(8.87) 


Now we have 

QlJp+,p.,M 


= 9} de 

•'P+,P_,0 


however, 
and similarly 

Together with 


2 5 (cos(6' - 9)) sin^(9 - 9) - \ cos{9 -9)\j MQ{p+xt + P-Xg ) 
+2D cos{9 - 9) sin(6' - 9) dg(^Mg{p+xt + P-X^)') 


S (cos(6' — 9)J Mg = 0, 

de(^Mg{p+x'l + P-X~g)) = {p+xt + P-X§) deMg. 
sm{9-9) 


( 8 . 88 ) 


dgMg{9) = -K 


cos‘^{9 — 9) 


Sign(cos(6* — 9y)Mg{9), 


we hnd (noting that kD = 1): 




cos(6* — 9) +2 


sin 2 (P-P) 


{p+xt + P-Xg ) Mg 


cos (6* — 

We hrst compute (the computation is obviously similar for T^). We have 


T+ = 


' cos(0—0)>O 


Ql\fp+,p.:e{0) d9 


= -n 


do 


P+,P_, 


' cos(0—0)>O 


cos(6* — 9) +2 


sin2(P - 9) 
{0-~9)\ 


cos 


{P+Xt + P-Xg ) M-g 


d9 


= -9\ 


■^P+,P_, 


cos(0 — 6*) I +2 


sin 2 (P-P) 
cos(6' — 9) 


{p+x§ + P-X-g ) M-g 


6/-0=+7r/2 


= 0 , 


6/-0=-7r/2 


because the contact of Mg to 0 at 9 — 9 = ±7r/2 is inhnite order and compensates for the 
blow up I cos(6* — 0)1 at the denominator. 
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We are left with the computation of T 2 . Using (3.33) and (8.88), we can write: 


To = 




= -e\ 




= -e\ 


''p+,p_, 


= -e\ 


•/p+,p_, 



t 

/ dg 


J — TV 

\ 

nTV 

r / 

/ dg 


/ —TV 

L\ 

nTV 

r / 

/ dg 


/ —TV 

L\ 


cos(6' — 9) \ +2 


sin 2 ( 0 - 0 ) 

W^\ 


cos 


{P+Xt + P-Xg ) M-g 


gg{9) d9 


cos 6* I + 2 
cos 6* I + 2 


sin 2 0 
I cos 6^1 
sinH 
I cos ^1 


{p+X^ + p-X )M 


g{9)d9 


M 


9i^) {p+X^ + P-X ) d9. 


We have 


dg 


sin^6*\ ,^1 ^ sin6* \ PsinP 

cos^ + 2- -— M = < Signfcos^) ( sin^ + 2 —— ) - k (-- + 

costly cos^9/ \cos9 


^ sin^ 9 
cos^ 9 


M. 


We note that this expression is changed into its opposite under the change 9 ^ 71 — 9. 
Since g M is also changed into its opposite under this transformation, it follows that the 
factors of p+ and p_ are equal, which leads to: 


To = -9\ 


''P+,P_, 


{sm9g)M + 2( 


cos^ 9 1m 

sin 9 


-K 


cos 9 /M 


; sin^ 9 

+ 2( —^ g , 

cos-^ 9 / M 


(P+ + P-)- 


Finally, inserting (8.87) into this expression leads to (5.58), which ends the proof. 
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